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We developed a novel universal eigenvalue analysis for 2D arbitrary nanostructures comprising dispersive 
and lossy materials. The complex dispersion relation (or complex Bloch band structure) of a metallic grating 
is rigorously calculated by the proposed algorithm with the finite-difference implementation. The abnormally 
large group velocity is observed at a plasmonic band edge with a large attenuation constant. Interestingly, we 
found the abnormal group velocity is caused by the leaky (radiation) loss not by metallic absorption (Ohmic) 
loss. The periodically modulated surface of the grating significantly modifies the original dispersion relation 
of the semi-infinite dielectric-metal structure and induces the extraordinarily large group velocity, which is 
different from the near-zero group velocity at photonic band edge. The work is fundamentally important to the 
design of plasmonic nanostructures. © 2015 Optical Society of America 
OCIS codes: 250.5403, 050.1755, 310.6628. 


For photonic crystals [1], forward and backward trav¬ 
eling waves constructively interfere with each other, re¬ 
sulting in standing waves. Due to Bragg scattering in pe¬ 
riodic photonic structures with modulated refractive in¬ 
dices, the degenerate standing waves split into two band 
edge modes with fields concentrated in different regions. 
Between band edge frequencies, a band gap occurs due 
to destructive interferences of traveling waves. Although 
it is nonintuitive that the phenomena could happen in 
metal structures, no physical law forbids their existences. 
Different from traveling waves in photonic crystals, for¬ 
ward and backward surface plasmon (SP) waves could 
also interfere with each other, leading to the formation 
of plasmonic band gap (PBG) and plasmonic band edge 
(PBE) [2]. They have broad applications in nanophoton¬ 
ics such as thin-film solar cells [3], lasers [4,5], surface- 
enhanced Raman scattering [6,7], directional nanoanten¬ 
nas [8], etc. 

The complex dispersion relation (or complex Bloch 
band structure) of metallic nanostructures plays a funda¬ 
mental role in understanding the PBE and PBG effects. 
In literatures, the dispersion relation has been approxi¬ 
mately generated by the frequency-angle diagram [6,9], 
which is not accurate for highly lossy metal materi¬ 
als. Because the eigenvalue (Bloch wave number) cor¬ 
responding to the incident angle will be complex given 
a real frequency. In this work, we developed a novel and 
rigorous eigenvalue analysis for 2D arbitrary nanostruc¬ 
tures with dispersive and lossy materials. Previously, the 
complex band structure can be calculated by converting 
the quadratic eigenvalue problem to linear one [10-12]. 
Here we borrow an idea from the quantum transport 
problem [13] to overcome the drawback. Eirst, the pro¬ 
posed method can handle universal dispersive and lossy 
media with experimentally-tabulated complex refractive 
indices. Second, the method can produce a linear eigen¬ 


value equation and thus saves computer resources in con¬ 
trast to the quadratic eigenvalue one. Third, it is easy 
to be implemented without decompositions of electro¬ 
magnetic fields. Einally, the proposed method can be 
employed to study the eigenvalue problem of arbitrarily 
shaped plasmonic nanostructures, which shows advan¬ 
tages over mode matching approaches [14-16] existing 
the numerical convergence issue. 

Using the proposed eigenvalue algorithm, we observed 
an abnormally large group velocity at the PBE of a 2D 
metallic grating with ID periodicity. The extraordinary 
group velocity distinguishes from the near-zero group ve¬ 
locity investigated both in photonic crystal structures 
[1] and in other periodic plasmonic structures with un¬ 
touched or separated unit elements [2,5]. 

Eig. 1(a) shows the schematic pattern of a squarely 
modulated metallic grating. SP waves will be excited at 
the crest and trough of the grating. The incident light 
is p polarized and vertically impinges on the grating. 
The complex refractive index of the metal (Ag) is ex¬ 
pressed by the Brendel-Bormann (B-B) model [17]. The 
mismatched momentum is provided by Eloquet modes 
of the periodic grating, which play an important role in 
inducing the PBE. Eig. 2(a) illustrates the optical ab¬ 
sorption of the metallic grating. The middle absorption 
peak is related to the PBG, while others are related to 
PBEs. The weaker (middle) absorption peak by the PBG 
is due to the coupling between top and bottom SP waves 
as shown in Eig. 2(c). When the thickness of the grat¬ 
ing T increases, the middle absorption peak disappears 
and thus an absorption dip is formed by the PBG (See 
Eig. 2(a)) where the reflectance has a high peak (See 
Eig. 3(a)). Eig. 3(b) depicts the phase distribution of Hz 
field at the PBG. We see a plasmonic standing wave pat¬ 
tern corresponding to the second-order Eloquet mode. 
While the planar semi-infinite air-Ag structure only sup- 
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Fig. 1. (a) The schematic pattern of a squarely modu¬ 
lated metallic grating with the geometric sizes of W = 
410 nm, P = 820 nm, H = 30 nm, and T = 50 nm. The 
red ellipses and yellow arrows denote the excited surface 
plasmons and incident light, respectively. (b,c) The nor¬ 
malized spatial overlap integrals between the near-held 
prohle of the grating and the m-th Floquet mode respec¬ 
tively at 410 nm (Peak 1 of Fig. 2(a)) and 480 nm (Peak 
3 of Fig. 2(a)). 
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Fig. 2. (a) The optical absorption of the metallic grating 
with a varying thickness T. (b,c,d) Hz held distributions 
of the metallic grating corresponding to the absorption 
peaks denoted by arrows (T = 50 nm). 



Fig. 3. (a) The rehectance of the metallic grating with 
a varying thickness T. (b) The phase distribution of Hz 
held at the PEG with respect to the Peak 2 in Fig. 2(a). 


To clarify the modihed dispersion relation, we theo¬ 
retically study the spatial overlap integral between the 
near-held prohle of the grating as shown in Figs. 2(b,d) 
and the m-th Floquet mode according to the rigorous 
coupled-wave analysis [18], i.e. 


T] = 

kx ' 



Hz{x, y) exp + jky{y - yc)) dxdy 


( 1 ) 

fco sin0 + 27rn/P, k^ + ky = kl, Im(fcy) < 0 (2) 


where P is the periodicity of the grating, and yc is the 
boundary between the grating crest and air. Moreover, 
is the wave number in free space, 6> = 0 for the 
vertical incidence case, and n = 0,±1,±2 are integers. 
From a simple calculation (Re(/co) > 27:jP)^ we know 
the zeroth- and hrst-order Floquet modes are radiative 
modes related to incident wave, specular rehection, etc. 
They are not responsible for the excitation of SPs. Fur¬ 
thermore, the overlap integral values become very small 
when the dihraction order |n| > 3 (See Figs. l(b,c)). 
Remarkably, the overlap integral value with respect to 
the third-order Floquet mode suddenly increases at 410 
nm as comparing Fig. 1(b) to Fig. 1(c). This also can 
be observed from the near-held distributions as depicted 
in Figs. 2(b,d), where the H-held is concentrated on the 
crest and trough of the grating, respectively. The inter¬ 
play of second- and third-order Floquet modes not only 
enables a coupling of SPs to photon energy but also per¬ 
turbs the original dispersion relation of SPs supported in 
the semi-inhnite air-Ag structure. A proper momentum 
matching condition can be revised as 


ksp + A/csp = k^ sin6> + 


(3) 


ports one plasmonic resonance peak that is signihcantly 
smaller than 400 nm, the two resonant peaks of PBEs at 
410 nm and 480 nm are obtained after introducing the 
grating structure. Intuitively, the original dispersion re¬ 
lation of the semi-inhnite air-Ag structure is modihed by 
the periodically modulated surface (refractive indices) of 
the metallic grating. 


where kgp = feoy , eo and are permittivities of 

free space and metal, is the perturbed momentum 

by the third-order Floquet mode, and n = 2 corresponds 
to the second-order Floquet mode. 

Besides the overlap integral method, the eigenvalue 
analysis is an efficient tool to characterize the PBE 
and PBG. Take ID periodicity for an example, the 
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wave equation can be discretized by using the finite- 
difference method and Bloch theorem (p^x -f- P) = 
(p{x)ex.p{-jkBP), i.e. 
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where p = [0i; (^ 2 ; ^ 4 ] is the eigenmode, Da = 

—2/A^ + [k{xyf‘^ T = 1/A^, and A is the spatial step. 


We define two new matrices H and Q as 
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(5) 

and use them to 

rewrite 

Eq. 

(4) as 








H/j) 





(6) 


Because Q containing columns filled with zeros cannot 
be inverted and H removing a diagonal element P 44 is 
singular, we should modify Eq. ( 6 ) to generate a well- 
conditioned linear eigenvalue problem 

fl-7\ 

(H-Qr'Q4'= ,I 

Mp = Xp. 

The matrix H — Q restores a standard finite-difference 
matrix of the wave equation and therefore can always be 
inverted. The above procedure can be easily extended to 
calculate the dispersion relation of the 2D metallic grat¬ 
ing with ID periodicity. The stretched-coordinate per¬ 
fectly matched layers [19] at the top and bottom bound¬ 
aries are adopted to absorb outgoing waves reflected by 
the metallic grating. At the left and right boundaries, 
the same Bloch theorem is employed as Eq. (4). It is 
worth mentioning the dispersion analysis of 2D gratings 
is an unbounded problem, which is fundamentally dif¬ 
ferent from 2D periodic crystal structures with bounded 
boundaries. 

The most commonly used method [1] for calculating 
dispersion curves is to choose the Bloch wave number 
ks beforehand, and the frequency is then computed as 
a function of the wave number, yielding the dispersion 
curve uj = uj{kB)- However, this traditional method is 
difficult to calculate the dispersion curve for the disper¬ 
sive material. In this paper, we fix the frequency uj and 
solve for the Bloch wave number as a function of fre¬ 
quency, i.e. kB = kB{oj). To benchmark the proposed 
method, we compute the band diagram of a periodic di¬ 
electric strip shielded with a perfect electric conductor 
(PEC) plate. Both methods agree with each other very 
well, which can be seen in Eig. 4. 



Eig. 4. The band structure of a periodic dielectric strip 
shielded with a perfect electric conductor plate. The di¬ 
electric constant and the side length of each strip is 4 
and P/2, respectively. 



Eig. 5. The band structure of the 2D metallic grating as 
a function of the real part of the Bloch wave number. 
A is the incident wavelength and the right-most verti¬ 
cal line denotes Re(/cB) = tt/P. The imaginary part of 
(physically real) metal permittivity by the B-B model 
increases or decreases by a factor of 2 respectively for 
simulating the high or low Ohmic loss. 

Next, we will investigate the complex Bloch band 
structure of a 2 D metallic grating as shown in Eigs. 5 
and 6 . At the PBE (denoted by the Peak I), the original 
dispersion relation of the semi-infinite air-Ag structure is 
modified after the introduction of the periodically mod¬ 
ulated refractive indices as illustrated in Eig. 5. More¬ 
over, we find an extraordinarily large group velocity at 
the PBE. Erom Eig. 6 , a large attenuation coefficient is 
found at the PBE corresponding to the large imaginary 
part of the Bloch wave number. To unveil the physical 
origin of the anomalous group velocity, we modify the 
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Fig. 6. The band structure of the 2D metallic grating 
as a function of the imaginary part of the Bloch wave 
number. A is the incident wavelength. The imaginary 
part of (physically real) metal permittivity by the B-B 
model increases or decreases by a factor of 2 respectively 
for simulating the high or low Ohmic loss. 

imaginary part of metal permittivity. Interestingly, the 
negative and normal group velocity are observed respec¬ 
tively for the low and high metallic absorption (Ohmic) 
loss (See Fig. 5). Hence, the leaky (radiation) loss leads 
to the extraordinary group velocity while the Ohmic loss 
smoothens the plasmonic band structure losing distinct 
features. The PBE, also called PBG edge or edge of PBG 
in literatures [8,20], is spectrally located far from the 
PBG in contrast to photonic band edge where the group 
velocity is almost zero. The interplay of second-order and 
third-order Floquet modes in the grating strongly per¬ 
turbs the original dispersion relation of the semi-infinite 
planar structure. The single modulated surface of the 
grating not only induces PBEs but also produces large 
leaky loss, which differs from plasmonic Bragg waveg¬ 
uides [14] where optical energy is confined between two 
modulated metal interfaces. 

The PBG at the peak 2 with a small attenuation co¬ 
efficient (See Eig. 6) is result from two counteracting 
SP waves with different group velocities (See Eig. 5). 
Gompared to PBEs with significant plasmonic hot spots 
at grating corners, the PBG has a weaker optical ab¬ 
sorption as depicted in Eig. 2(a). In other words, the 
PBG suppresses the Ohmic loss. The PBG is located at 
Re{kB) = 0 that is the edge of irreducible Brillouin zone 
similar to the photonic band gap. It should be noted 
that the zero PBG is caused by the symmetric property 
of metallic gratings and the PBG could open up in dual¬ 
interface gratings [20]. 

The work proposed a novel universal eigenvalue solu¬ 
tion to investigate the dispersion relation of a metallic 
grating. The extraordinary group velocity at the PBE 
is confirmed by our eigenvalue solver. The interaction 
between Eloquet modes supported by the periodically 
modulated surface induces the strong leaky loss and thus 


the abnormally large group velocity. Both the devel¬ 
oped numerical algorithm and new findings with detailed 
physical explanations are useful for designing plasmonic 
nanostructures. 

The authors acknowledge the support of the grants 
(Nos. 712010, 711609, and 711511) from the Research 
Grant Gouncil of the Hong Kong, the grant (No. 
10401466) from the University Grant Gouncil (UGG) 
of the University of Hong Kong (HKU), and the Seed 
Project Eunding of HKU (No. 201211159147). This 
project is also supported by the UGG of Hong Kong (No. 
AoE/P-04/08), by The National Natural Science Eoun- 
dation of Ghina (No. 61201122), and in part by a Hong 
Kong UGG Special Equipment Grant (SEG HKU09). 


References 

1. J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. 
D. Meade, Photonic Crystals: Molding the Flow of Light 
(Princeton University Press, 2008). 

2. D. de Ceglia, M. A. Vincenti, M. Scalora, N. Akozbek, 
and M. J. Bloemer, AIP Adv. 1, 032151 (2011). 

3. X. H. Li, W. E. I. Sha, W. C. H. Choy, D. D. S. Fung, 
and F. X. Xie, J. Phys. Chem. C 116, 7200-7206 (2012). 

4. T. Okamoto, F. H’Dhili, and S. Kawata, Appl. Phys. 
Lett. 85, 3968-3970 (2004). 

5. W. Zhou, M. Dridi, J. Y. Suh, C. H. Kim, D. T. Co, 
M. R. Wasielewski, G. C. Schatz, and T. W. Odom, Nat 
Nano 8, 506-511 (2013). 

6. A. Kocabas, G. Ertas, S. S. Senlik, and A. Aydinli, Opt. 
Express 16, 12469-12477 (2008). 

7. N. Mattiucci, G. D’Aguanno, H. O. Everitt, J. V. Fore¬ 
man, J. M. Callahan, M. C. Buncick, and M. J. Bloemer, 
Opt. Express 20, 1868-1877 (2012). 

8. M. Toma, K. Toma, P. Adam, J. Homola, W. Knoll, and 

J. Dostalek, Opt. Express 20, 14042-14053 (2012). 

9. W. L. Barnes, T. W. Preist, S. C. Kitson, and J. R. 
Sambles, Phys. Rev. B 54, 6227-6244 (1996). 

10. A. Raman, and S. H. Fan, Phys. Rev. Lett. 104, 087401 

( 2010 ). 

11. C. Fietz, Y. Urzhumov, and G. Shvets, Opt. Express 19, 
19027-19041 (2011). 

12. M. Luo and Q. H. Liu, J. Opt. Soc. Am. A. 27, 1878- 
1884 (2010). 

13. M. Luisier, A. Schenk, W. Fichtner, and G. Klimeck, 
Phys. Rev. B 74, 205323 (2006). 

14. C. Li, Y.-S. Zhou, and H.-Yu. Wang, Opt. Express 20, 
7726-7740 (2012). 

15. T. Okamoto and S. Kawata, Opt. Express 20, 5168-5177 
( 2012 ). 

16. R.-L. Chern, C. C. Chang, and C. C. Chang, Phys. Rev. 
E 73, 036605 (2006). 

17. A. D. Rakic, A. B. Djurisic, J. M. Elazar, and M. L. 
Majewski, Appl. Optics 37, 5271-5283 (1998). 

18. M. G. Moharam, E. B. Grann, D. A. Pommet, and T. 

K. Gaylord, J. Opt. Soc. Am. A. 12, 1068-1076 (1995). 

19. W. C. Chew, and W. H. Weedon, Microw. Opt. Technol. 
Lett. 7, 599-604 (1994). 

20. T. Okamoto, J. Simonen, and S. Kawata, Phys. Rev. B 
77, 115425 (2008). 


4 









